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Abstract 

With one Am^ of the appropriate order of magnitude to solve the atmo- 
spheric neutrino problem, we study the resulting three-generation vacuum- 
oscillation fit to the solar neutrino flux. An explanation of the atmospheric 
neutrino composition in terms of pure — > Vr oscillations is easily compat- 
ible with the well-known two-flavour oscillation solution for solar neutrinos. 
The allowed parameter region in the sin^(20i2)-Am^ plane changes little 
with increasing values of the mixing element C/gS , provided this is less than 
about 0.4. We find that the threefold maximal mixing is disfavoured. 
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1 Introduction 



While neutrino-oscillation experiments tend to suffer from systematic errors that 
are difficult to estimate, there are at least two distinct neutrino problems that 
deserve serious consideration: the deficit of solar neutrinos and the atmospheric 
neutrino composition. The former has been around for many years. It started 
with the pioneering Homestake detector and later gained significance with the 
water-scattering and gallium experiments [1]. 

The latter anomaly is an inconsistency between the number of detected neutri- 
nos of the muon family and those of the electron family, assumed to be generated 
by cosmic-ray hadrons impinging on the atmosphere. The ratio of ratios 
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where iV^ is the number of /i-like events and Ne the number of e-like events, 
should be, according to theory, close to 1, but is instead around ~ 0.6 [2]. 

If one analyses the solar- neutrino deficit (a summary of the chlorine [3], gal- 
lium [4, 5] and water-scattering experiments [6, 7] is given in Table 1) in terms of 
two-flavour mixing, the vacuum-oscillation solution leads to a set of disconnected 
regions in the two-dimensional parameter space spanned by sin^(2^) and Am^, 
according to the Ug survival probability 

P..^..(t) = l-sin^(2^)sin^(^), (2) 

where 9 is the mixing angle, Am^ the difference of the squared masses, t the time 
since creation and E the neutrino energy. The favoured region is sin^(20) ~ 0.8- 
0.9, with Am^ ~ 6 x 10^^^ eV^, but neighbouring regions that differ from Am^ 
by n X SmF' can yield comparable values of x^, provided n is a small integer, and 
Sm^ is of the order of 1.4 x 10~^^ eV^ (cf. Fig. la)^. The scale (5m^ is given by 
the characteristics of the detectors, the neutrino energy (of the order of a MeV) 
and the distance R from the Sun to the Earth: Srri^ ~ AnEhc/R. 

The probability expression (2) and Fig. la show the 'just so' character of 
vacuum oscillations, the mass-squared difference has to be fine-tuned in relation 
to the Sun-Earth distance if this model shall solve the solar- neutrino problem [9]. 
(The allowed region in the upper part of Fig. la, where Am^ > 2 x lO"^'^ eV^, is 
a new feature of the most recent data [10].) 

With three flavours of neutrinos, one may consider the general three-flavour 
mixing, in which a state, which is of pure flavour a {a = e, /i, r) at t = 0, evolves 
as 

ry{t)^Y.Uaje-'^^'iy„ (3) 

j 

^Thc MSW [8] interpretation of the data leads to values of Am^ that are several orders of 
magnitude higher. 
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where the summation index j runs over mass states. We write the unitary mixing 
matrix as 

Uel Ue2 Ue3 

U^i U,2 . (4) 

Url Ut2 Ut-s 

The full parameter space then consists of three angles (we ignore CP- violating 
effects) and two independent squared-mass differences (the energy can be approx- 
imated as Ej ~ p + m'j /2E). In particular, the survival probability of an electron 
neutrino takes the form 
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where Amf^- = mf — m'j. Unitarity and data on other phenomena, in particular 
on the atmospheric neutrino flux, lead to constraints on these parameters. 

We denote the mass states ui, and 1/3, according to their intrinsic mass 
ratio: their respective mass eigenvalues will be sorted as mi < m2 < 777.3. For the 
mixing matrix we use the form advocated for quarks by the Particle Data Group 

[11]': 



U 



U^r{923)Uer{9l3)UeM 

C12 Ci3 S12 Ci3 Si3 

~Ci2 Sl3 S23 ~ ^12 C23 "512 Si3 S23 -|- C12 C23 Ci3 S23 

— C12 Sl3 C23 + S12 S23 "512 Si3 C23 — C12 S23 Ci3 C23 



(6) 



where C12 = cos ^12, C13 = cos ^13, etc., and the mixing angles ^12, ^13 and ^23 are 
in the interval [0, 7r/2]. 



2 Three-parameter accommodation of the solar 
neutrino flux 

A three-flavour analysis of the solar-neutrino oscillations will restrict the allowed 
values of the elements in the mixing matrix. Because of unitarity, the solar 
electron-neutrino survival probability can always be expressed in terms of only 
two mixing angles, here taken to be 612 and ^13. We shall assume that the dom- 
inant mixing (represented in Fig. la) is a mixing between the first two families, 
and ask (i) 'How will the goodness of fit depend on C/es' and (ii) 'How will the 
allowed regions of Fig. la change as one increases the strength of this mixing?' 

Before we address these questions, we want to introduce some further restric- 
tions on the neutrino parameters: we reserve one /\m? for the vacuum-oscillation 

^ There are nine structurally different ways to express the mixing matrix in terms of mixing 
angles. A discussion of the parametrization of flavour mixing in the quark sector is presented 
in [12], where also CP violation is included. 
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solution of the solar-neutrino problem, Am|i ~ 10~^^ eV^. The suppression of 
atmospheric muon neutrinos relative to electron neutrinos suggests that also 
oscillates into some other flavour of neutrinos, but (since the ratio t/E relevant 
to atmospheric neutrinos is orders of magnitudes smaller than that for solar neu- 
trinos) with a mass-squared difference that is many orders of magnitude larger 
than is required to explain the solar-neutrino deficit. It is therefore natural to 
assume that other squared-mass differences are much larger: 

Am^^ ~ 10-1^ eV^ < Am^2 ^ (7) 

such that when Eq. (5) is applied to the case of solar neutrinos, the last two terms 
can be represented by their average values: 

—^\--sm\2e,,). (8) 

Here, the last term represents a non-zero component of the u-^ mass state in the 
electron- neutrino state, and vanishes with |C/e3p. For ^13 = we get the familiar 
two- family mixing of Eq. (2). 

While a simultaneous flt to both solar- and atmospheric-neutrino fluxes also 
would involve the angle 6'23 as well as two mass-squared differences, under the 
assumption (7) the above formula (8) provides a three-parameter accommodation 
of the solar fiux. 

We show in Fig. lb the results of two-parameter fits to the solar-neutrino data. 
(The high region in Am^, cf. Fig. la, is not displayed here.) In our calculations 
we use the solar-neutrino fiux as given by the 1998 Standard Solar Model (SSM) 
by Bahcall, Basu and Pinsonneault [10], and for the detection rates the values 
given in Table 1. The second mixing angle, ^13, is held fixed at different values, 
and the fit is performed in the variables sin^(2^i2) and Am|^. It is seen that 
regions that are separated in the two-family case, can be smoothed out to one 
big region when the more general case of three fiavours is considered. However, 
there is a slight worsening of the quality of the fit as the coupling to the third 
state is increased. For sin ^13 = 0.0, 0.2, and 0.4, the values of Xmin ^'^^ 3.6, 4.2 
and 5.6, respectively. For two degrees of freedom, 4 experiments (SAGE and 
GALLEX are combined) minus 2 parameters, the probability of getting larger 
than those values by chance, are 17, 12 and 6%, respectively. None of these is 
particularly good, but they are not excluded. 

Since there is some uncertainty associated with the important boron fiux 
[10], we have checked how the fit changes under 10% rescalings of the fiux. In 
Table 2 we show the resulting x^j^ and goodness of fit for three different fractions 
of the SSM boron flux (keeping the other components of the flux fixed), for two 
different values of Ue^. For the lowest of these fractions, the goodness of fit is 
appreciably better than for the standard ^B-fiux. It is also seen that, for the 
lowest of our tabulated boron fiux values, the fit deteriorates less when sin ^13 is 
non-zero. 
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A different representation, for several values of sin ^13, is given in Fig. 2, 
including also the higher- Am^ region. As sin ^13 increases, the allowed regions 
get broader, but also the fit gets poorer. 

It has been argued that the data suggest 'maximal mixing' [13]. In our vari- 
ables, maximal mixing corresponds to sin^(26'i2) = 1 and sin ^13 = l/\/3 ~ 0.58. 
We find that the fit gets rather poor for such parameters, with Xmin = 11- 

3 Long-baseline experiments 

The future long-baseline experiments [14] will provide new restrictions on these 
mixing elements. For example, if we assume that a long-baseline experiment is 
such that oscillations other than those involving 

AM^ = Amgi ~ Am32 (9) 

can be neglected, then the condition (7) leads to the following transition proba- 
bilities: 

Pit.M - sin^ 023 sin^(2^i3) sin^ , (10) 

Pit.At)^ cos' e^s sin' (20,,) sin' , (11) 

and the survival probability: 

P!;^.M) = 1 - ^-^(2^13) sin^ (^) • (12) 

Note that the coefficient in Eq. (10) can be written as 4|C/e3P|t^//3|^- The 
transition can thus be given a simple physical interpretation, valid for small values 
of the mixing angle 6*13. When ^m^^t/E <^ 1, the transition can be interpreted 
as an 'indirect' transition, involving the admixtures of as given by ^^13 and 
^23- This is different from the case of Am\^t/ E ~ 1, where the 'direct' transition 
~^ is effective via the states vi and 1^2. 

The atmospheric-neutrino data are usually explained in terms of Ue or 

^fj. oscillations with a large amplitude [15, 16]. When applying Eq. (8) to 

explain the solar-neutrino deficit, the transition z/^ is preferable for the 

atmospheric neutrinos, because it allows for a small or vanishing value of 6'i3, 
which gives the best fit of the observed solar neutrino rate. An early three- 
fiavour analysis (not including data from Super-Kamiokande) gave for the best 
fit a value for \Ue?}\' (or sin^^^i3) in the range 0.06-0.14 [17]. More recent studies 
(including Super-Kamiokande data) prefer values for sin ^13 of 0.1 [15] and < 0.1 
[16]. As we have seen, these low values of sin 6*13 also lead to acceptable fits to 
the solar-neutrino data, with x^j^^ = 3.6-5.6. 

In the long-basehne experiments [14, 18], neutrinos from either an accelerator 
or a reactor will be counted after travelling a 'long' distance, from 1 km (for 
CHOOZ) to 730 km (for the MINOS experiment). In fact, the recent data from 
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the CHOOZ experiment, where the survival probabihty Pg is measured, does 
not cxchidc any of the cases studied in Fig. 2, provided AM^ < 9 x 10~^ eV^. 
Indeed, no evidence for oscillations was found for AAP > 2 x 10'^ eV^ and 
sin^(26'i3) > 0.18. This means that sin6'i3 has to be either small (< 0.22) or quite 
big (> 0.98) (or AM^ < 2 x lO^^ eV^). With the large-angle solution, Eq (8) 
gives a poor fit to the solar-neutrino data, this alternative should therefore be 
disregarded. The remaining low-angle solution is also preferred in the analysis of 
the atmospheric-neutrino data [15]. In this case (with AM2 > 2 X 10-2 eV^) the 
lower two boxes of Fig. 2 would be excluded. The Bugey reactor experiment [19], 
which is less sensitive, excludes AM^ down to about 3 x 10"^ eV^ for sin ^^13 > 
0.16, consistent with CHOOZ. 

A measurement of a non-zero value for Pj, ^^.^ will imply a non-zero value for 
^^13. Such information, which would be a major result from the planned long- 
baseline experiments, combined with Eq. (8), can help us to further restrict the 
allowed parameter space relevant to the solar neutrinos. 

We note that, for <-> oscillations, the K2K experiment [20] will have 
a sensitivity of sm'^{29) > 0.1 and the MINOS experiment could probe values 
as low as sm'^{29) = 0.01. In our parametrization, this means that they would 
be sensitive to values of sin ^13 as low as 0.16 and 0.05. If this mixing element 
should turn out to be significantly larger, the vacuum-oscillation solution to the 
solar-neutrino problem would have serious difficulties. 

4 Concluding remarks 

The data from the solar-neutrino detectors can very well be reconciled with the 
atmospheric-neutrino data. If we assume that the solution of the latter anomaly 
is the — > transition (which is induced by the direct mixing 6*23 and gives the 
best fit), then the vacuum-oscillation solution for the atmospheric and the solar- 
neutrino problems can be regarded as two separate two-generation oscillations. 
Even if the transition z/^ — ^ z/g, which for Am\^L/E ^ 1 is induced by the two- 
step process of v^^r and v^Vt mixing (with angles ^^13 and ^^23, respectively) should 
turn out to take place, the rate has to be quite large before it can severely alter 
the fit of the vacuum-oscillation solution to the solar-neutrino problem. 

When the fgi^r mixing angle ^13 is increased from zero, the allowed ranges 
of Arn^^ merge to a wider one, but beyond sin^^i3 ~ 0.5, which is well within 
the sensitivity that will be reached by the K2K experiment, the fit deteriorates 
rapidly. 
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lixi)c>iiiiic>iit 


Counting rate 


ooi\ r r i nl 
ooi\i [KJJ 


Lountcxl/boM 


HoTTiPstakp f3l 


2 54 ± 14 ± 14 


7 7+ln 
' ■ ' —1.0 


330 ± 026 


GALLEX [4] 


69.7 ±6.7111 


1291^ 


U.O+U_o 063 


SAGE [5] 


79+12+5 
' ^-10-7 


1291^ 


U.OOO_o.095 


Kamiokande [6] 


2.80 ±0.19 ±0.33 


c 1 c:+0.98 
^•-■-^-0.72 


0.544 ±0.074 


Super-Kamiokande [7] 


2.42 ± 0.06l[j:^? 


c 1 c:+0.98 
"-'•-'-^-0.72 


0.4701H?? 



Table 1: Solar neutrino fluxes. Measured and SSM- predicted event rates are for 
the radio-chemical detectors measured in units of SNU (1 event per second per 
10^^ target atoms). For the scattering detectors (the last two entries) the fluxes 
are in units of 10^ cm^^ s^^. In the ratio of the observed to the predicted values, 
only experimental uncertainties are included. 





sin ^13 = 0.0 


sin ^^13 = 0.3 


/b = 0.9 /b = 1.0 /b = 1.1 


/b = 0.9 /b = 1.0 /b = 1.1 


A mm 


2.0 3.6 5.8 

38 17 5.5 


2.9 4.8 7.2 

23 9 3 



Table 2: Best flts for different fractions, /b, of the SSM flux, for two different 
values of sin ^13 = Ues- The lowest row gives the goodness of fit. 
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0.5t , , , — = ^^^^ 

0.7 0.8 0.9 1 

sin2(2^i2) 




Figure 1: Allowed parameter space within 95% CL. The upper figure illustrates 
the two-generation case, the cross marks the position of Xmin- The lower figure 
shows (under the assumption leading to Eq. (8)) the dependence on the mixing 
element t/es = sin ^13, which increases in steps of 0.1 from 0.1 to 0.4. The regions 
merge as the third mass state i^s gets more involved in the mixing. 
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0.5 0.75 

sin2(2^i2) 



2.5 



CD d. 



1.5 



0.5 




sin6'i3 = 0.4 




0.5 0.75 1 

sin2(2^i2) 



Figure 2: Allowed parameter space within 95% CL. for AM^ ^ Am^, for dif- 
ferent values of t/es (= sin6'i3). In the two upper boxes, the allowed regions for 
the two-generation oscillation (Fig. la) are also shown, as shaded areas. Note 
different horizontal scales. The crosses show the best-fit point. 
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